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ABSTRACT 

Late X-ray flares observed in X-ray afterglows of gamma-ray bursts (GRBs) suggest late 
central engine activities at a few minuets to hours after the burst. A few unambiguously con- 
firmed cases of supernova associations with nearby long GRBs imply that an accompanying 
supernova-like component might be a common feature in all long GRB events. These motivate 
us to study the interactions of a late jet, responsible for a x-ray flare, with various components 
in a stellar explosion, responsible for a GRB. These components include a supernova shell-like 
ejecta, and a cocoon that was produced when the main jet producing the GRB itself was prop- 
agating through the progenitor star. We find that the interaction between the late jet and the 
supernova ejecta may produce a luminous (up to 10 49 erg s _1 ) thermal X-ray transient last- 
ing for ~ 10 s The interaction between the late jet and the cocoon produces synchrotron-self 
absorbed non-thermal emission, with the observed peak X-ray flux density from 0.001 /iJy to 
1 mJy at 1 keV and a peak optical flux density from 0.01 /iJy to 0. 1 Jy (for a redshift z = 2). 
The light curve due to the late jet - cocoon interaction has very small pulse-width-to-time 
ratio, At ft ps 0.01 — 0.5, where t is the pulse peak time since the burst trigger. Identifying 
these features in current and future observations would open a new frontier in the study of 
GRB progenitor stars. 
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1 INTRODUCTION 

Long duration gamma-ray bursts (GRB) - those lasting more than 2 
s - are thought to be produced by a relativistic outflow (or jet) with 
a kinetic energy ~ 10 51 — 10 52 erg (beaming effect corrected) when 
a massive star collapses at the end of its nuclear burning life cycle 
(see Piran 1999, 2005; Meszaros 2002 for reviews). The massive 
star origin of GRBs are supported by two different observations: 
(i) GRBs are found to be in actively star forming galaxies (e.g., 
Christensen et al. 2004, Castro Ceron et al. 2006) or in star (es- 
pecially massive ones) forming regions of the host galaxies (e.g., 
Paczyriski 1998; Bloom, Kulkarni & Djorgovski 2002; Fruchter et 
al. 2006); (ii) for a subset of nearly a dozen of GRBs, X-ray-rich 
GRBs and X-ray flashes, Supernova (SNa) features - both tempo- 
rally and spacially associated with the bursts (Woosley & Bloom 
2006 for a review) - were detected. For four of those: GRB 980425 
(e.g., Galama et al. 1998), 030329 (e.g., Hjorth et al. 2003), 02121 1 
(Delia Valle et al. 2003) and 031203 (e.g., Malesani et al., 2004), 
the physically associated SNe were not only photometrically but 
also spectroscopically confirmed. The others of the subset show a 
late-time (~ 10 days) SNa-like "bump" in the optical afterglows, 
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with a simultaneous strong color evolution, e.g., in GRB 980326 
(Bloom et al. 1999) and 011121 (Bloom et al. 2002), consistent 
with the hypothesis of an underlying SNa. 

The Swift satellite has recently unveiled a "canonical" behav- 
ior pattern in about two-thirds of GRBs' early X-ray afterglows: a 
rapid decline phase lasting for ~ 10 2 s is followed by a shallow 
decay phase lasting ~ 10 3 — 10 4 s, then by a "normal" power-law 
decay phase and finally by a possible jet break (Nousek et al. 2006; 
O'Brien et al. 2006). In addition, X-ray flares are found in about 
50% of all Swift bursts; they have been discovered in all of the 
above phases (Burrows et al. 2005, 2007; Chincarini et al. 2007). 
Even long before Swift, a late X-ray flare was detected by Bep- 
poSAX for GRB 970508 (Piro et al. 1998). 

In this work we investigate a scenario in which a late jet - 
responsible for producing the X-ray flares and possibly the shal- 
low decay phase - interacts with other components of a long-GRB 
stellar explosion. These components include a SNa ejecta, if a SNa 
explosion is accompanying the GRB event, and a cocoon created 
by the passage of the main GRB jet through the star (Ramirez-Ruiz 
et al. 2002, Matzner 2003, Zhang et al. 2004). 

In Section 2 we argue for the existence of the late jet and mul- 
tiple components of a GRB explosion, and provide physical moti- 
vations of this work. We investigate the interactions of the late jet 
with the SNa ejecta in Section 3 and with the cocoon in Section 4 
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and calculate their associated emissions. The predicted emissions 
and their detection prospects are confronted with current observa- 
tional data in Section 5. The summary and implications are given 
in Section 6. 



2 A MULTI-COMPONENT GRB STELLAR EXPLOSION 

Two recently discovered GRB features point to the emergence of a 
late outflow (jet) after the main 7-ray producing outflow has died. 
The first is the X-ray flares observed at a few x 10 2 — 10 3 s (as 
late as 10 4 — 10 5 s in some cases) after the prompt burst, with 
a fluence typically about one tenth of the fluence of the prompt 
7-rays; in one case, GRB 050502B, this ratio is ~ 1 (Falcone et 
al. 2007). Lazzati, Perna & Begelman (2008) found that the mean 
flux of a flare declines with its occurrence time as ~ t~ 1,5 . These 
flares are characterized by a large flux increase and by a very steep 
rise and decay. Typical increase of the flux ranges from a factor of 
order unity to 10 and in some rare cases even a few hundred. The 
decay after the flare peak is as steep as oc t~ 4 , much steeper than 
the underlying afterglow decay (oc t^ 1 ). The pulse width to the 
peak time ratio At/t is much smaller than unity, typically ~ 0.3 
(Burrows et al. 2005, 2007; Chincarini et al. 2007). 

The late flares cannot be due to external-origin mechanisms 
such as, a density clump in the circumburst medium, or the energy 
injection into the afterglow blastwave by the trailing slower shells, 
because the decay slope after the rebrightening from these mecha- 
nisms always follows the standard afterglow model, and At/t ~ 1 
is always expected (Nakar, Piran & Granot 2003; Nakar & Gra- 
nd 2007; Lazzati & Perna 2007). They are also unlikely to arise 
from late collisions between two slow moving shells ejected at 
same time as the main 7-ray producing shells, since the resultant 
internal shock is too weak to give rise to the significant emission 
observed in flares (Zhang 2006). The most likely possibility for X- 
ray flares is the late activity of the central engine (e.g., Fan & Wei 
2005). Such a late activity was proposed already in 1998 as an al- 
ternative origin of GRB afterglow (Katz & Piran 1998, Katz, Piran 
& Sari 1998). 

The scenario involving the late engine activity can easily sat- 
isfy the constraint that At/t <C 1. Also in this scenario the late 
flare is physically separate from the "background" afterglow, so the 
large amplitude increase of the flux superposed on the decaying 
afterglow can be naturally explained. 

The second feature favoring the existence of a late jet is the 
shallowly decaying component in the overall X-ray light curve. 
This component starts at a few xlO 2 — 10 3 s and ends at 10 4 — 10 s s 
during which the flux decays slower than expected for standard de- 
celerating blast-wave afterglow. One straightforward interpretation 
is that the shallowness is due to an long-lasting energy injection 
from a late outflow that catches up with the decelerated main out- 
flow (Cohen & Piran 1999; Zhang & Meszaros 2001; Yu & Dai 
2007), though there are other possible interpretations for this fea- 
ture (e.g., Granot & Kumar 2006; Granot, Konigl & Piran 2006; 
Fan & Piran 2006, Panaitescu et al. 2006; Kobayashi & Zhang 
2007; Uhm & Beloborodov 2007; Genet et al. 2007). 

In those cases where we don't see a SNa spectroscopically or 
photometrically in the optical band, we can still hope to explore 
the existence and properties of the SNa accompanying the GRB by 
looking at the interaction of a late jet with the SNa remnant and the 
emission from it. It is likely that every long-duration GRB has a 
SNa accompanying it but some extrinsic and/or intrinsic bias might 
have hindered the optical detection of the SNa component(Woosley 



& Bloom 2006). A late jet provides a chance to test this picture in 
the cases where the ordinary SNa features are not easy to detect. 

Given the evidences for the existence of a late jet and the 
physical association of SNe and GRBs, one expects a number of 
different interactions between the following four components in a 
GRB event (see Fig. 1 for an illustration): (i) A highly relativistic 
(r ~ 300 — 1000) narrow jet along the rotation axis with an open- 
ing angle ~ 0. 1 radian that produces the main GRB event; (ii) A 
nearly spherically symmetric SNa ejecta moving with speed ~ 10 4 
km s ; (iii) A cocoon fireball created by the passage of the main 
GRB jet through the star; (iv) A late relativistic jet responsible for 
the late X-ray flares after the end of the main GRB. See Woosley & 
Heger (2006) for a detailed version of the multi-component GRB 
scenario. Also see Wheeler et al. (2000) for a similar but more de- 
tailed model which relates various energetic phenomena such as 
SNa, GRB and magnetar in a single stellar explosion. 

In this work we investigate the interactions of a late jet which 
might have a similar Lorentz factor and opening angle to those of 
the main GRB jet, with the expanding SNa ejecta and the cocoon. 
We assume that the SNa ejecta and the main GRB jet are launched 
from the central source at about the same time . Wheeler et al. 
(2002) described a similar stellar collapsing scenario where a de- 
layed relativistic jet from the eventually formed black hole catches 
up and collides with an earlier proto-pulsar toroidal field generated 
mildly relativistic jet as an origin of the 7-ray burst. Ghisellini et al. 
(2007) firstly considered the collision between a jetted fireball from 
an intermittent GRB central engine and a stationary cocoon as an 
alternative to the standard internal shock scenario, but with higher 
efficiency. Here we are not attempting to explain how a relativistic 
jet is formed and a 7-ray burst is produced from the relativistic jet 
or how an accompanying SNa is generated, rather we are trying to 
constrain the physical properties of the late jet, SNa ejecta and the 
cocoon by calculating the emissions from their interactions, and to 
verify the general picture of GRB-SNa connection. 

We will use the following fiducial values for a variety of model 
parameters. For the late jet, a total energy Ej ~ 10 51 erg, opening 
angle 9j ~ 0. 1 radian, LF Vj ~ 100, a delay respective to the 
launching of the main GRB jet, cocoon and SNa ejecta (in our pic- 
ture the three are launched at more or less the same time) tF ~ 10 2 
s, and a duration tdur ~ 10 2 s; for the SNa ejecta, we use an 
isotropic-equivalent mass Msn ~ 10 Mq and speed Vsjv ~ 10 9 
cm s -1 ; for the cocoon, we use an energy E c ~ 10 51 erg, a ter- 
minal LF r c ~ 10 and an opening angle 6 C ~ 0.6 radian. Nev- 
ertheless an appropriate range of numerical values is assigned to 
each parameter in the real calculation (e.g., Ej could be two orders 
of magnitude larger or smaller than the fiducial one, cocoon speed 
could be sub-relativistic and t F could be as large as 10 4 — 10 5 s). 

The late jet will run into the SNa ejecta first (at a distance ~ 
10 1 cm; see Fig. 2a), and then catches up and run into the cocoon - 
if it successfully crosses the ejecta - at a larger distance (~ 10 12 — 
10 14 cm; see Fig. 2b). 



1 The observational constraint is that GRB and SNa occur within ~ 1 day 
(e.g., Woosley & Bloom 2006). SNa explosion theories estimate that the 
SNa shock breaks out at a few tens of seconds after the core rebounce (e.g., 
Janka et al. 2007 and references therein); this time scale is small compared 
to the delay of the late jet (t F > 10 2 s). 
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Figure 1. Schematic illustration of multiple components in a long GRB and an accompanying supernova (SNa). An initial, highly relativistic jet, as shown 
in the right end of the illustration, is responsible for the prompt GRB. A cocoon, that was inflated by the initial GRB jet as the jet was punching through the 
envelope of the progenitor star, has broken out of the stellar surface at the same time as the main jet, and has accelerated to mildly relativistic speed. A nearly 
spherically symmetric sub-relativistic ejecta is responsible for the SNa. A late jet responsible for the late X-ray flares in GRBs is launched from the central 
source and will catch up with the above components. The "late" cocoon is produced when the late jet crosses the SNa ejecta. The distances are not to scale. 



3 LATE JET - SNA EJECTA INTERACTION 

Since the late jet is highly relativistic and the SNa ejecta is sub- 
relativistic, the late jet will catch up with the SNa ejecta in a time 
roughly equal to the distance of the SNa ejecta from the explo- 
sion center (where the jet emerged), vsn, divided by light speed 
c. If the late jet is launched with a delay of tp, then vsn ~ 
VsNtF ~ 10 VsN,9tF,2 cm (hereafter we use the convention 
X n = X/10 n unless specifically notified). At this time the ra- 
dial width of the ejecta is about the same size of the initial stellar 
envelope, Agjv ~ r* ~ 10 11 cm. The ejecta has hardly moved 
from the initial position of the progenitor stellar envelope for a 
jet delay tp ~ 100 s. The particle density in the SNa ejecta is 

n S N ^ 10 24 Vg^ gtp^ 2 M S N,i^sN,u cm ~ 3 > where M S n is in 
units of the solar mass, and it is extremely optically thick. If there 
is any emission from the interaction between the late jet and the SN 
ejecta, that emission should be thermal. 

3.1 The cavity in the SNa ejecta 

We recall that the main GRB jet has already traversed the star and 
left a cavity in the polar region before the late jet comes. Since the 
material enclosing the "wall" of the cavity was heated up by the 
passage of the main jet, it has tendency to refill the cavity. To find 
out whether the cavity in the SNa ejecta has been filled up before 
the late-jet encounter, we estimate the time it takes for filling up the 
cavity. 

When inside the star, the cocoon material has a relativistic 
temperature, i.e., the local sound speed c s , c = c/y/3. Thus the 
cocoon material may fill the cavity on a time scale of r t 9j/c StC ~ 
1 s, much shorter than the onset of the late jet ti?. However the 
cocoon will also break out and flow away from the ejecta at that 
same speed in ~ 10 s, leaving behind a somewhat evacuated polar 
region. The filling of the polar cavity by the rest of the SN ejecta 
is uncertain. Assuming a temperature T ~ 10 s K for the SNa- 



shocked ejecta material, it has c s ~ 1.2 x 10 8 T g 1,/2 cm s 1 , so the 
filling time would be r*8j /c s (T) ~ 10 a r* ; iiTg s, comparable 
to tf. Considering that the ejecta local temperature possibly de- 
creases from inner parts to outer parts and the transverse size of the 
cavity gets bigger outward, it is likely that when the late jet hits the 
ejecta the cavity is partly filled - the inner part is filled but the outer 
part is not. However, for a very late jet, e.g., tp ~ 10 3 s, the cavity 
is surely filled before the late jet comes up; and then the question 
that arises is whether the jet is powerful enough to cross the new 
refilled cavity. This issue will be discussed elsewhere. 

3.2 Late jet - SNa ejecta crossing 

The late jet comoving particle density is rij = 
Lj/inryri^mp) = 7 x 10 16 L^rr^^er^ 
cm -3 , where Lj is the late jet luminosity. Because usn /nj ^> 1, 
the late jet will be decelerated to a non-relativistic speed after it 
first hits the ejecta. It will undergo the same process as the GRB 
main jet did when propagating through the progenitor star. A 
cocoon forms in the SNa ejecta and makes a way for the jet head 
by pushing the ejecta material sideways. To distinguish it from the 
cocoon associated with the main GRB jet, let us call this cocoon 
associated with the late jet as the "late cocoon" (see in Figs. 1 and 
2a). 

When the jet is moving inside the ejecta, the jet head has been 
slowed down to be at a sub-relativistic speed (Ramirez-Ruiz et al. 
2002; Matzner 2003): 

r h y /2 

= 0.8 x W^^Or^Ms^A 1 /^ cms" 1 . (1) 

Note that if the SNa ejecta width is constant then Vh has no de- 
pendence on tsn and therefore on tp; but if the SNa ejecta is uni- 
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Figure 2. (a): Schematic illustration of the late jet - SNa ejecta interaction, 
at a slightly later time than shown in Fig. 1. The breakout of the "late" 
cocoon produces a short thermal emission lasting for ~ 10 s. (b): Schematic 
illustration of the late jet - cocoon interaction, at a somewhat later time than 
in (a), when the late jet has completely emerged out of the SNa ejecta and 
is colliding with the cocoon created by the main GRB jet as it made its way 
through the progenitor star. A pair of shocks are going through the cocoon 
and the late jet, respectively. Note that the "late" cocoon has expanded and 
cooled down at this time. 



formly distributed from the centre to the radius vsn, i.e., Asn 

/2 

: . SNa exph 
rsN, is the most probable 



1/2 

rsN, then increases with tF as oc i F ' . SNa explosion Simula 



tions show the latter case, i.e., Asn 
one (e.g., Tanaka et al. 2009). 

A constraint on the late jet property can be derived from the 
requirement that the duration tdur of the jet must be larger than 
the time that the jet spends to cross the SNa ejecta. Thus tdur > 



1/2 

Asn/vh implies L/^tdur^O^ 



3.3 Thermal emission from the late cocoon break out 

The late cocoon - formed by the interaction of late jet with SNa 
ejecta 2 - will break out of the star immediately following the break- 



2 The late cocoon should not be confused with the cocoon that was formed 
and left behind by the main GRB jet. The latter will be discussed in Section 
4. 



out of head of the late jet (Fig. 2a). The luminosity of the thermal 
emission from the late cocoon breakout can be estimated from its 
temperature and transverse size. Prior to its breakout, the late co- 
coon has a pressure p c and its leading head moves with the jet head 
at the same speed in the radial direction and expands transversely 
into the SNa ejecta with a speed v± . The ram pressure balance at 
the lateral interface between the late cocoon and the ejecta material 
gives 



Pc = Psnv ± , 



(2) 



where psN is the ejecta mass density. 

The late cocoon contains an energy E c that is approxi- 
mately equal to the jet luminosity Lj times the shell crossing time 
Asn/vh, and it is radiation pressure dominated. The volume of the 
late cocoon just before the breakout is V c = TvAsNr^_/3, where 
r± = AsNV±/vh is the transverse size. So the pressure is 

E c LjVf, 
PC= Wc = 



ttA 2 



(3) 



Combining with Eq. (2) gives 

1/2 

Wh 1 

ttA 2 



_ ( L]PSNVh\ 



(4) 



and the thermal temperature is 



Tth 



V a 



1/4 



= 1.2 x 10 B L'gO^. 



3/16a-l/8 „, 1/16 -1/4 



-1 Ivl SN,l' SN^l^SN, 11 Jrv : 



(5) 



where Eq. (1) is used. For a late jet with ~ 10 s, the typical 
thermal photon energy should be a few keV. Tth becomes smaller 
for larger tF', for instance, when the SN ejecta width Asn is ~ 
rsN and, if Lj oc tp 1,5 (Lazzati et al. 2008), we have T t h oc tp°' g . 

Let us estimate the luminosity of the thermal emission from 
the late cocoon breakout. The late cocoon transverse expansion 
speed is 



PSN 



1/2 



U.o X 1U L , 49 tf • , M SJV1 r SJV11 A 



-1/8 



(6) 



The cocoon transverse size is r± = Asnv±/vh = 0.8 x 



10" £ 3 - ^^-^^Sff/Swi^,!! Cm ' N ° tiCe th£lt ^ iS al " 

most ~ rsN for fiducial parameter values. However, the trans- 
verse size of the visible emitting region of the cocoon just at the 
breakout should be smaller than r±. This is because when the jet 
head moves near the outer surface of the ejecta, it probably accel- 
erates due to the rapid density drop of the stellar materia there; 
thus the lateral expansion of the cocoon's leading head immedi- 
ately following the jet head should be suppressed. Therefore at that 
time, the overall cocoon might be in an "hourglass" shape, as il- 
lustrated in Ramirez-Ruiz et al. (2002) and in our Fig. 2a, rather 
than a conic shape. At the breakout, the hot cocoon material es- 
capes and accelerates radially from a "nozzle" which has a trans- 
verse size on the same order of the jet's. Thus, the transverse size of 
the visible emitting cocoon right at the breakout can be estimated 



3/8 



as r t h,± = rsN0j = 10 



-iVsN,9tp,2 cm, while r± should be 



the transverse size of cocoon at its widest cross section. 
The luminosity of the black body radiation is 

L t h = <?Tt h Trrt h:± 

* r 1r >48 7-3/4 a 3/2 wl/4 \ — 5/4 

= 4.5 x 10 L^e^Ms'^rsN.ix^sNM er S s 



(7) 
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Using the fact that Lj oc t F , one can find this thermal luminosity 



decreases with tF as oc t„ ' if A< 



rsN, or as oc t. 



if 



Asjv ~ constant. Thus, the radiation efficiency increases as oc t F 
for Asjv ~ rsjv and as oc t F 4 for Asjv ~ constant, respectively, 
which is simply because the emitting transverse area increases with 
t F . 

This thermal transient will last for a time comparable to the 
time it takes for the bulk of the late cocoon to escape the ejecta. Af- 
ter that, the luminosity drops quickly due to adiabatic cooling. The 
late cocoon's outflow speed is ~ c s = c/v3, so the escape time is 
At e3C m Agjv/c s = 6 Asjv, 11 s, Note that Lth is large because 
all the jet power was deposited in the cocoon over a relatively long 
time (ejecta crossing time ~ 10 2 s), and this energy is efficiently 
radiated away (via black body radiation) within a relatively short 
time, ~ 6 s. 

An even stronger (L ~ 10 49 erg s _1 ) thermal pulse is associ- 
ated with the main GRB jet break out. However, this thermal tran- 
sient is short (~ 10 s), with spectral peak at X-rays, and it happens 
during the initial stage of the main "burst" 7-ray emission when 
the X-Ray Telescope (XRT) is not pointing towards the burst. Since 
this emission can be estimated to be one to two orders of magni- 
tude dimmer than the emission of the burst itself, it is difficult to 
observe. On the other hand, the thermal transient due to the late co- 
coon breakout that we consider here arises later and at the time that 
the XRT is already pointing towards the burst. Hence this transient 
should be easier to detect. 

The fact that we don't observe such thermal transient, preced- 
ing flares during the afterglow is somewhat puzzling. If the late jet 
is fairly early (say, tp *S 10 s), the polar cavity in the SNa ejecta 
left by the main GRB jet was probably only partially filled, i.e., 
only the inner part of the cavity is filled (see Sec. 3.1). In that case, 
the luminosity of the thermal emission would be much smaller than 
what we estimated above. The situation is more problematic for 
tF ^ 10 li . We discuss further the observational prospects of this 
transient in Sec. 5.1. 



4 LATE JET - COCOON INTERACTION 

We now turn to the interaction of the late jet, after it has success- 
fully crossed the SNa ejecta, with the cocoon that was formed by 
the main GRB jet (see Fig. 2b). The cocoon breaks out from the 
star at the same time when the main jet breaks out. Then it accel- 
erates to a mildly relativistic speed. The delay of the late jet with 
respect to the cocoon breakout is tF and the late jet catches up with 
the adiabatically cooled cocoon at n — ctp/{Pj — j3 c ), where /3j 
(Tj) and /3 C (T c ) are the speed (LFs) of the late jet and the co- 
coon, respectively. When the cocoon speed is mildly relativistic, 
n » 6X 10 14 t F 2r 2 1 cm, which is much further than the late jet - 
SNa ejecta interaction site. For a sub-relativistic cocoon, n is close 
to but still outside the jet - SNa ejecta interaction region. 



4.1 Cocoon geometry and dynamics 

At the breakout, the cocoon has an energy E c , and an energy-to- 
mass ratio ?7 C . The cocoon opening angle 6 C is determined by its 
transverse expansion speed ~ c/v3, thus C ~ I/a/3 = 0.6. 
During the early stages of expansion, the cocoon's radial width, 
A c , is approximately the width of the stellar envelope r*. Later 
on, due to the radial expansion of a relativistically moving gas, 
A c asymptotically approaches r/(2r 2 ) in the lab frame; this hap- 



r c (f) ~ 8 c r/r„ when r < r s , and T c f» r\ c when r ^ r s , where 
r s = T] c r*/6 C is the saturation radius (Paczynski 1986; Good- 
man 1986; Shemi & Piran 1990; Piran, Shemi & Narayan 1993; 
Meszaros, Laguna & Rees 1993). The evolution of the cocoon's 
comoving volume, Vc'(r) = ir0 2 r 2 A c r c (r), is described by 



V' 



Tr^r 3 , 
n8 2 r 2 ri c r,, 



for r < r s 

for r s < r < r u 

for r > r w . 



(8) 



The pressure evolution of the cocoon follows the adiabatic 
expansion law: p c oc V c ~ 7 . Initially the radiation pressure domi- 
nates, so 7 = 4/3. When the cocoon's optical depth to Thomson 
scattering decreases to below unity, the photons decouple from the 
matter and the radiation pressure drops exponentially; then the gas 
pressure takes over the dominance of the pressure with 7 = 5/3. 
The transition happens at radius rt — [oTE c /(ri c mpC 2 n9 2 )] 1112 = 
3.5 x W 14 El^ricl^ 2 cm, where ot is the Thomson scattering 
cross section. The initial pressure at the breakout is given by 



= 5 x 10 17 Ecfiir'u dyn cm -2 . 



" ' 3Vc'(r.) 
The evolution of the cocoon pressure is thus given by 

k r ) > 

'9,^/3 /r, \8/3 



Pc(r) 
Pc,0 



(2r, c e c ) 



4/3 



for r < r s , 
for r 3 < r <r w , 
for r > r w . 



(9) 



(10) 



The comoving density of the cocoon is assumed to be homo- 
geneous; the same is for the late jet. The width of the late jet is 
determined by its duration and the radial expansion, so Aj (r) = 
ctdur + T "/(2r 2 ). The jet comoving density is 



rij(r) 



E, 



4nr 2 m p c 3 r 2 (t dur + -£-)' 



(11) 



pens when r ^ r u 



The cocoon's LF can be described as 



4.2 Cavity in the cocoon 

There was also initially a cavity in the cocoon left by the passage 
of the main jet. Here we estimate how quickly the cavity would be 
filled. The filling up process starts when the main jet dies off. For 
a typical duration of GRB t gr b ~ 10 s, the cocoon has moved to a 
radius r c ~ ct gr b — 3 x 10 11 t flr 6,i cm. Since r c < r s , the cocoon 
gas is still relativistic, the sound speed is c a ~ c/y3, and the time 
required for cavity to close is ~ r c 6j/c 3 — 0.6 t 9r b, i#j,-i swhich 
is <C tF - Thus the cavity is securely filled when the late jet reaches 
the cocoon. 

Now let us consider a possibility that there is a continuous 
low-level central engine activity (a low-power jet) following the 
end of the main GRB jet and preceding the late jet. This low-power 
jet and the late jet more or less are the same phenomenon, only with 
different energy fluxes. The low-power jet might have too small 
radiation luminosity to have an observational imprint, however it 
might still be dynamically important for keeping open the polar 
cavity in the SNa ejecta and the cocoon. 

To keep the cavity open without significant energy dissipation 
from the jet and cocoon, the transverse pressure, i.e., the thermal 
pressure, of the low-power jet should be greater than or equal to 
the cocoon's pressure at r c . Though the low-power jet started with 
a high thermal pressure at a distance ro ~ 10 2 km from the ex- 
plosion centre where it was launched, at r c its thermal pressure has 
dropped to be much smaller than the cocoon pressure because it 
has adiabatically expanded by a much larger factor than the cocoon 



© 0000 RAS, MNRAS 000, 000-000 



6 Shen, Kumar & Piran 



does. Thus, when this low-power jet entered the cavity left by the 
main GRB jet, it will be squashed rapidly by the gas pressure in the 
cocoon. 

The ram pressure of the low power jet, Pj, r am(j c ) — 
Lj,iow I '(^nr^c) where Lj.i ow is the jet luminosity, can help bore 
a cavity through the cocoon under suitable conditions. For a long- 
lasting, continuous jet, the event that takes place after it has been 
squashed is as follows. If Pj. ra m{r c ) > p c (r c ), then the resid- 
ual, incoming jet will punch a new channel through the cocoon. 
In doing this, the jet is heated up by the reverse shock, so that 
the jet which is moving inside the channel could have an en- 
hanced pj (thermal pressure) that is comparable top c , and can keep 
the channel open. If we assume the low-power jet was launched 
at the end of the GRB, i.e., t gr b ~ 10 s, making the equality 
Pi,ram = Pc(r c ) gives Lj tlow = 2 x 10 48 erg s _1 , with other 
parameters at fiducial values. The same equality at later times gives 
the time dependence Lj > i ow {t) oc t~ 2 . So a decaying luminosity 
profile Lj > i ow (t) = 2 x 10 4S t~[ 2 erg s^ 1 for the low-power jet is 
required to keep the cavity open, i.e., the minimum required total 
energy for a low-power jet to keep the cavity open is a few x 10 49 
erg. Note that at this minimal jet luminosity, the process of keep- 
ing the cavity open is not smooth and a fraction of jet energy is 
dissipated and it should have some radiation associated with it. 

In conclusion, the polar cavity left by the main GRB jet in 
the cocoon can fill up quickly - before the late jet reaches the co- 
coon. The presence of a continuous low-power jet can keep this 
cavity open, provided that the jet has a minimal total energy of a 
few x 10 49 erg. In the following calculation we consider the case 
when the cavity is filled up. Clearly if the cavity is empty the late 
jet cocoon interaction is trivial. 

4.3 The jet - cocoon interaction phases 

The dynamical process of the late jet - cocoon interaction can be 
decomposed into three phases in the following sequence 3 : (i) The 
Collision Phase takes place when the jet runs into the cocoon with 
a forward shock (FS) propagating into the cocoon and a reverse 
shock (RS) propagating into the jet (see Fig. 2b). (ii) The Penetra- 
tion Phase begins when either the RS crosses the entire jet or the FS 
crosses the entire cocoon, whichever comes first. In the first case 
the shocked fluid (RS-shocked jet and FS-shocked cocoon fluid) 
will decelerate, , after RS crosses the jet, when more and more co- 
coon material are swept by the FS. In the latter case the shocked 
fluid will be accelerated, after FS crosses the entire cocoon, by the 
remaining unshocked jet ejecta, and a new particle population will 
be accelerated at the RS. (iii) The Expansion Phase begins when 
both the FS and RS have run through their courses, and the entire 
shocked fluid expands adiabatically. 



4.3.1 The collision phase 

In this phase, the FS propagates into the cocoon and the RS prop- 
agates into the jet. For simplicity we approximate the interaction 
using a planar geometry. The entire jet / cocoon system can be di- 
vided into several zones. Outside the FS (RS) front is the unshocked 
cocoon (jet), and these are taken to be cold plasma, i.e., e = n and 
p = 0, where e, p and n are the fluid energy density (including the 

3 The three-phase decomposition treatment follows Dermer (2008) who 
studied the emission due to the external shocks between a GRB jet and a 
stationary circumburst dense cloud. 



rest mass energy), pressure and particle number density, respec- 
tively, all measured in its comoving frame. In between the FS front 
and the RS front are the shocked cocoon fluid and the shocked jet 
fluid. These shocked fluids move with the same LF (r s ) and have 
the same thermal pressure; they are separated by a contact discon- 
tinuity (CD) plane. 

The fluid properties across a shock front are governed by the 
mass, momentum and energy conservation laws (e.g., Landau & 
Lifshitz 1959; Blandford & McKee 1976). Across the shock the 
fluid particle density increases by a factor of (7F + — 1), 

and e = Tn, where T is the shocked fluid LF measured in the 
unshocked fluid comoving frame; 7 is given by: p — (7— l)(e— n). 
We use m p — c — 1 to simplify the formulae. 

At the CD, the pressure in the shocked jet material equals the 
pressure in the shocked cocoon material, i.e., 



(if 



i)(f s 



(7* 



\ IRS 

s-i)(r sc - 



+ 1 



V 1FS 



1 



(12) 



where the subscript "RS" refers to the reverse-shocked fluid and 
"FS" the forward-shocked one, and f s j and r sc are the shocked 
fluid LFs measured in the unshocked jet and cocoon comoving 
frames, respectively. 7 lies between 4/3 and 5/3, and can be written 
in terms of f as 7 = (4T + l)/3f (Kumar & Granot 2003). Then 
Eq. (12) simplifies to 

{f% - l)n 3 = (f L - l)n c . (13) 

Since L SJ = TjT^l-pjP,) and f sc = L s r c (l - p s p c ), we find 
from the last equation that 



T s = r, 



(14) 



J (a + 1 + 2^/ET jc y/ 2 ' 

where a = nj/n c is the density ratio, Fj C is the unshocked jet 
LF measured in the unshocked cocoon rest frame. Note that this 
expression for T s is valid for both sub-relativistic and relativistic 
shocks, and for 4/3 ^ 7 ^ 5/3. 

The shock (RS or FS) front moves with a LF different from 
the LF of the shocked fluid. The shock front LF as measured in 
the unshocked fluid comoving frame - which we denote as Trsj 
for the RS and as Tfs,c for the FS - is given by the solution to the 
shock-jump conditions as a function of 7 and T (Eq. 5 of Blandford 
& McKee 1976). Using the expression of 7 in terms of T, we find 



RS,j 



41- - 1 



VST 



and Tfs. 



4H 



Vsric + 1 



(15) 



These expressions are valid for both sub-relativistic and relativistic 
shocks. 

The pair of shocks exist until one of the two shocks, either the 
RS or FS, has traversed through the unshocked fluid. From that time 
on, the interaction will move to the next dynamic phase (Penetra- 
tion). To determine which shock (RS or FS) crossing occurs first, 
we calculate two radii, r^s an d r^s, where r^s is the distance of 
the system when the RS crosses the rear end of the jet, and tfs is 
when the FS crosses the front end of the cocoon, pretending that 
the pair of shocks had existed all the way to the larger of the two 
radii. 

At trs or rps, the total distance that the shock has traveled 
through the unshocked fluid is equal to the radial width of the jet 
or the cocoon at that radius. Thus the two radii can be obtained by 
solving the equations 
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Figure 3. Left panels: the comparison between the RS crossing radius r^g and the FS crossing radius rpg in the late jet - cocoon interaction for a set of model 
parameter space. Right panels: the comparison of the two crossing radii with 2r^, where ri is the radius at which the late jet catches up with the cocoon and 
the interaction begins. T c and (3 C are the cocoon's LF and dimensionless speed, respectively. Ej and E c are the kinetic energies of the late jet and the cocoon, 
respectively. Both crossing radii are numerically calculated from Eq. (16). Other parameter values are: E c = 10 50 erg, 9j = 0.1, Q c = 0.6, t F = 300 s, 



Aj (rns) = Wj ~^ RS) dr and 

MrFs)^!" 3 iPFS ~/ c) dr, (16) 

where we have used dt = dr/ fij ~ dr//3 c . 

If rns < rFS, the RS crosses the jet before the FS crosses 
the cocoon, and vice versa. We calculate trs and tfs for different 
parameters. The results are shown in Fig. 3. We find that rus < 
vfs, i-e., the RS crossing occurs first, for most of the parameter 
space; rns > rFS can only happen when the cocoon bulk motion 
is sub-relativistic (T c f3 c < 1) and the energy carried by the late jet 
is much larger than that of cocoon (Ej / E c 3> 1). We will use this 
result later (Sub-sections 4.4 - 4.5) to simplify the calculation of 
the light curve by assuming that RS crossing always occurs before 
FS crossing. 

4.3.2 The penetration phase 

After the RS crosses the jet, the FS continues to pass through the 
cocoon. We consider the entire shocked fluid, both the old and the 
newly shocked, moving together with the same LF. In the rest frame 
of the unshocked cocoon, the LF of the shocked fluid is determined 
by the equation for the deceleration of a relativistic blast wave in the 
adiabatic limit as the blast wave sweeps up the stationary ambient 
medium: 

r'jv') = — — r ^- A ( 17) 



(Bottcher & Dermer 2000) where the prime sign denotes the un- 
shocked cocoon rest frame, T' s A is the shocked fluid LF at the end 
of the collision phase, m c (x') is the isotropic equivalent swept-up 
mass at the distance x' that the blast wave has traveled. 

A significant deceleration of FS occurs after a point where 
F's ' A m c(x' d )c 2 = E'j iso . Before this point, the FS is coasting into 
the unshocked cocoon at roughly the same speed it had prior to the 
RS crossing. After this point, the FS decelerates as it sweeps more 
and more cocoon material (similar to the external shock scenario 
for the GRB afterglows). If the FS crossing is earlier than the RS 
crossing, the process is similar except that it is the RS that continues 
to travel through the unshocked jet. 

4.3.3 The expansion phase 

After the FS eventually crosses the entire cocoon, the shocked fluid 
expands outward with a LF determined by eq. (17) but with m(x') 
replaced by the isotropic equivalent total mass of the cocoon. The 
relativistic electrons cool via radiation and expansion. The radia- 
tion from an adiabatically expanding relativistic shell, when it is 
optically thin to the Thomson scattering, is discussed in Barniol 
Duran & Kumar (2009). We will address the optical thick case in 
Sec. 4.5.1. 

4.4 The emission from the late jet - cocoon interaction 

We calculate the emission from the late jet - cocoon interaction, 
and estimate the flux densities at the optical and X-ray bands using 
standard shock synchrotron emission (e.g., Sari et al. 1998), and 
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Figure 4. The observed spectrum from the late jet - cocoon interaction at the 



peak of the light curve for Ej 



10 50 i 



3, tdur/tp — 0.3 



and 2 = 2. Other parameter values are same as in Fig. 3. Various lines 
represent contributions from different emission regions and spectral com- 
ponents: red dashed - RS synchrotron; red dotted - RS SSC; blue dashed - 
FS synchrotron; blue dotted - FS SSC; solid - the sum. For these parameter 
values, the order of the frequencies is v c < v op t < v a < < vx ; fo P t 
is in the synchrotron-self-absorption optically thick regime, and the spectral 
peak is around UV band. The SSC emission contribution is important only 
for 7-ray band and above but is negligible at both optical and X-ray bands. 



taking into account the synchrotron self Compton (SSC) radiation. 
Just behind the shock front, a fraction of the bulk kinetic energy of 
the fluid upstream, e e , is transferred to the electrons, and another 
fraction, e B , goes to the magnetic field. All the swept-up electrons 
are shock-heated into a power-law energy distribution with a spec- 
tral index p. The minimum LF of shock heated electrons is 



7i 



P 



(r-i), 



(18) 



where F is the shocked fluid LF measured in the unshocked fluid 
frame. The magnetic field energy density downstream of the shock 
front is given by 

B' 2 

U' B = — = AV{Y -l)e B n m v c 2 , (19) 
air 

where B' is the comoving frame field strength and no is the particle 
number density of the unshocked cocoon or jet. The synchrotron 
characteristic frequency that corresponds to ji is 

* = 9 (\1 V (20) 

27rm e c(l + z) 

where e is the electron charge. 

The electron cooling LF, 7 C , is determined by radiative cooling 
due to synchrotron and SSC radiations: 

j c m e c 2 = ^a T cj 2 c U' B (l + Y)t', (21) 

where ctt is the electron's Thomson cross section, t is the elapsing 
time in the shocked fluid rest frame, and Y is the Compton param- 
eter defined as the ratio of the SSC to synchrotron luminosities; 7 C 
is obtained by numerically solving Eq. (21) (e.g., McMahon, Ku- 
mar & Piran 2006). Electrons with LF > j c will cool to 7 C in time 
t'\ the cooling effect is negligible for electrons with LF < j c . The 
synchrotron cooling frequency is 

r s eB'j 2 c 



2nm e c(l + z) 



(22) 



The self-absorption frequency for the synchrotron electrons, 
v' a , measured in the shocked fluid comoving frame is calculated by 
(Sari & Piran 1999; Li & Song 2004;, McMahon et al. 2006; Shen 
& Zhang 2009) 



max(7 m ,7 a ) x 2m e v' 2 = F' v , , 



(23) 



where F' v , is the flux density at v' a radiated away from the surface 
of the shocked region. 

The emergent synchrotron spectrum of the shock-heated elec- 
trons can be approximated as a piece-wise power law function. The 



peak of the /„ -spectrum is at v n 
density at the peak is 



N e T a m e c 2 a T B' 



i'. -ma ,r 



= min(^, v c ) and the flux 



(24) 



4ttD 2 3e(l + z) ' 

where N e is the isotropic equivalent total number of shock-heated 
electrons, D is the luminosity distance. We also calculate the ob- 
served flux density due to SSC by (Rybicki & Lightman 1979) 
3 

} %C {v) = -OT&S 



(25) 



where Ss is the line-of-sight width of the emitting source, v s and 
f svn {v s ) are the synchrotron frequency and flux density (in the 
observer frame), respectively, and n e (-y) is the number of shocked 
electrons per unit volume per unit interval of 7; n e (7) oc 7~ p for 
7 > 7i [the modification of n e (-y) due to the radiative cooling is 
included in our calculations]. The function F(x) = 2x lnx + x + 
1 — 2x 2 for < x < 1 and is otherwise. Using the expression for 
optical depth r e due to Thomson scattering, the SSC flux density 
can be written as 



vt y 7i 7 2 ^47 2 u 3 



(26) 



Fig. 4 depicts the spectrum observed at the time when the RS 
crosses the jet. The order of the characteristic frequencies in this 
example is v c < u opt < v a < Vi < vx , and u a is in the UV band. 
The SSC contribution is negligible at optical band and it becomes 
important for photon energy 1 keV. 



4.5 Light curves 

The light curve from the the late jet - cocoon interaction is mainly 
determined by the evolution of f Vj7nax , Vi, v c and v a - The opti- 
cal depth of the shocked fluid region may alter the final light curve 
shape, which will be addressed later in this sub-section (§4.5.1). 
We follow the treatment of Yu & Dai (2009) and calculate the light 
curve. We define T exp as the time when the shocked fluid has trav- 
eled a distance of 2r, , where is the interaction radius (here and 
in the following, times denoted with the capitalized letter "T" are 
the observer's times and T — is the time when the interaction 
begins). Thus, T exp — Ti/{2Y 2 c); before T exp , the increase of the 
radius can be neglected and B' and F 3 are considered to be con- 
stant; after T exp , the attenuation of the density and B' due to the 
radius increase must be taken into account. 

We also define the shock-crossing time T cro = 
min(Tns ,Tfs), where Trs and Tfs, calculated in Eq. (16), 
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Figure 5. Contours of r e — the shocked region optical depth to Thomson 
scattering at the peak of the flux from the late jet - cocoon interaction. Top: 
for tp = 10 2 s; Bottom: for tp = 10 3 s. For even larger tp (say ~ 10 4 
s), we find r e < f for all parameter space. Other parameter values are same 
as in Fig. 3 except t^ ur /tp = 0.3. 



o 



o 



LO 



LO 

o 



LO 

o 



T F g < 1 





10 < T FS < 100 



> 100 



10 ° 0.01 0.1 1 

E./E 

J c 




T FS < 1 




10 



io < 

^3 



0.01 



HHP* 

0.1 1 
E./E 



100 



10 



100 



Figure 6. Contours of 77? s — the optical depth for the unshocked cocoon 
that is still in front of the FS at the peak of the flux from the late jet - 
cocoon interaction. Top: for tp = 10 2 s; Bottom: for tp = 10 3 s. For even 
larger tp (say ~ 10 4 s), we find rps < 1 for all parameter space. Other 
parameter values are same as in Fig. 3 except t^ ur /tp = 0.3. 



are the crossing times for the reverse shock and the forward 
shock, respectively (Fig. 3 shows Trs < Tfs for most of the 
model parameter space). Before T cro , the radial spreading of the 
shocked region is suppressed by the existence of two shocks, 
thus the volume of the shocked region V' oc r 2 and the internal 
energy density e' oc V'^ 1 oc r -2 ; the total number of shock 
heated particles increases linearly with time. After T cro , the radial 
expansion has to be considered and the shocked region experiences 
adiabatic cooling. During this phase, V' oc r" (where s = 2 ~ 3) 
and the internal energy density e! oc V'^ 4 ^ 3 oc r _4s ^ 3 . 

Therefore the evolution of B' and f v , m ax are as follows: (i) 
For T cro < 



B' 



T u , T < T CXF 

T -2s/3 rp > j, 



(27) 



T, T < T oxp ; 

a,x oc ^ T , 2exp "C T <C Tcro 5 
r -2s/3 ; r>Tcro . 



(30) 



And the temporal dependences of the characteristic frequen- 
cies are given by 



T 
T~ 



T < T ra 



-2s/3 



T < T ox 



for T cro < T cxp , and 



T° T <r T 

± 5 ' ' C X p 5 

T , T cxp ^ T *v T cro , 
r -2s/3 ; T>Tcro; 



(3D 
(32) 

(33) 



T, T<T cro ; 

/: i/,max OC ^ .£ , Tcro "n T < ?exp j 

T- 2s/3 , T>T oxp ; 

(ii) For T cro > T cxp , 

T°, T < T cxp ; 
B' oc <( T-\ T oxp < T < T CIO ; 
T~ 2s/3 , T>T CIO ; 



T <T c: 



(28) 



(29) 



(34) 



for T cro > 

T cxp . 

Note that the evolution of v c given above ignores SSC cooling. 
SSC cooling is included in our numerical calculation as described 
by Eq. (26). We find that the resultant v c scalings are not very dif- 
ferent from those given above. 

To calculate the light curve, we first calculate B', /„,„„, Vi, 
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Figure 7. The comparison between the observed in situ flux density at T p 
and the photospheric flux density, at the optical (top) and X-ray (bottom) 
bands, respectively. The photospheric flux is meaningful only in cases of 
(r e + rps) > 1 at Tp\ those cases where (r e + rps) < 1 are marked 
with fl h = 0. In cases of rps > 1, the emission produced at T p is unde- 
tectable, and the light curve is dominated by photosphere emission arriving 
at a later time; these cases are marked with fv{T p ) = 0. When r e > 1 and 
T FS < 1. me light curve is dominated by either the in situ emission pro- 
duced at T p and diminished by a factor of T e or the photospheric emission, 
whichever is larger, so is the observed peak flux. 



v c and v a at the expected peak time T p = min(T cro , T exp ), then 
we use the temporal evolution of B' and f v . m ax to get the observed 
flux density at other times. The SSC contribution to the flux is in- 
cluded in the light curve calculation. In addition, the optical thick- 
ness of the late jet - cocoon system to Thomson scattering could 
alter the light curve shape because it could delay the escape of the 
photons from the system; we will consider this next, and then dis- 
cuss the light curve results. 



4.5.1 Optically thick cocoon 

When n is small - either because of a low LF of the cocoon 
or a small tp - the cocoon could be optically thick to Thomson 
scattering. For instance, the optical depth of the entire cocoon at 
n is estimated to be « 2oTE c /[A-Kr1 F c m p c 2 (l — cos(9 c )] ~ 
0.1 E c ^\T~\t F 2 2 . Figs. 5 and 6 depict the calculated r e and rps — 
the optical depths of the shocked region and the unshocked cocoon, 
respectively — at the expected light curve peak time T p . 

When r e 1, the photons are subject to numerous scattering 
(diffusion) before escaping the plasma. The emergent flux is spread 
over the diffusion time scale 
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Figure 8. Top: temporal evolutions of B' , i/j, v c , v a and f VyT nax during the 
late jet - cocoon interaction for same model parameter values as in Fig. 4. 
Vi, v c and v a are in units of eV, whereas B' and f v , m ax are normalized by 
their maximum values, respectively. Bottom: observed light curves in opti- 
cal and X-ray bands. The observer's time T is normalized by the expansion 



time scale T e 



: n/{2TgC) and T = is when the interaction begins. 



AT d 



r p h 

2r§ c ' 



(35) 



which is the delay between the actually observed time of a photon 
and the time it would have been observed in the absence of scatter- 
ing, where r p h is the photosphere radius. 

When r e ^> 1, we consider all the photons emitted during the 
time up to T p as a photon gas co-expanding with the baryon gas; 
the expansion of the system is governed by the radiation pressure 
rather than the gas pressure, and the scattering between photons and 
electrons is nearly elastic. This treatment is different from the one 
adopted by Pe'er, Meszaros & Rees (2006) who consider the case 
where the gas pressure dominates over the radiation pressure and 
photons Compton scatter off thermal electrons. 

The equation of state for the photon gas is {hu) i oc V -4 / 3 , 
where hv is the characteristic photon energy and V oc r 2 A c is 
the volume of the system (in the photon-baryon gas co-expanding 
phase, the LF of the system is constant). The width of the system 
A c is oc r for the thin shell case and is constant for the thick shell 
case. Since the system is shock-compressed after the collision, we 
believe the thin shell case is a more likely possibility to consider 
than the thick shell case. Thus the temporal evolution of the photon 
energy is hv oc T _1 , where T is the observer's time. The shape of 
the spectrum at the photosphere is unchanged from that at the time 

We first calculate the flux density at T p neglecting all optical- 
thick effects. If r e > 1 and tps < 1, the promptly observed flux 
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Figure 9. Top: same as Fig. 8 but for Ej = 10 X E c = 10 51 erg and 
r c = 5. Bottom: observed light curves. Both the optical and X-ray light 
curves show a flat part at T < T exp because both frequencies are in the 
fast cooling regime and the decrease of u c oc T~ 2 just compensates for the 
increase of fv^max <x T. The optical peak at T > T exp is due to the fact 
that w c increases steeply and v a drops below v op t ■ A rise does not show up 
in the X-ray light curve at T > T exp because i>x is always above all of 
v a , v; and v„. 



in situ at T p is 1/ r e of that calculated when r e is neglected. Then 
we also calculate the emergent flux at the photosphere. The real 
peak flux of the light curve is either the observed flux in situ at T p 
or the photospheric flux at later time, whichever is larger, and the 
light curve would be dominated by that larger component. In case 
the observed flux in situ at T p is stronger than the photospheric 
flux, we calculate elaborately the light curve shape following T p 
in the way described in Sec 4.5 and numerically estimate its peak 
time and pulse width. In case the observed light curve is dominated 
by the flux at the photosphere, its peak time is the photospheric 
time T p h ~ r p h/(2F 2 c) and the pulse width is also ~ T p h- If 
i~fs > 1, then the unshocked cocoon blocks the light produced in 
the shocked region from reaching the observer, thus the emission at 
the photosphere is what we actually see only and it will completely 
dominate the light curve. We show in Fig. 7 a comparison between 
the observed flux in situ at T p and the photospheric flux for the 
considered model parameter space. 

4.6 Results 

Figs. 8-10 depict the light curves and the temporal evolutions 
of emission properties for different values of model parameters. 
Initially the light curve remains constant up to the shock crossing 
time or the expansion time, whichever is smaller. This is because 
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Figure 10. Top: same as Fig. 8 but for Ej = 10 2 X E c = lir* erg 
and r c /3 c = 0.5. Bottom: observed light curve in X-rays. The optical light 
curve is not shown, because for this set of parameter values the cocoon is 
extremely optically thick (r e 2> 1 and rp$ = 0; see Fig. 5 and 6), and 
we find the optical light curve is dominated by the photosphere emission, 
whose numerical light curve shape has to be calculated differently from 
that in the optically thin case. However in the X-ray band, the prompt non- 
thermal flux, after diminished by the optical thick effect, is still brighter than 
the later photospheric flux. Thus we use the diminished prompt non-thermal 
flux to represent the observed X-ray light curve. 



v c < v pt, and therefore the increase of f v ,max is compensated 
by the decrease of v c . Later, the X-ray flux decays due to the adia- 
batic expansion, while the optical flux continues to rise as long as 
v op t < v a - Therefore generally the optical pulse peaks later and is 
wider than the X-ray pulse. 

The peak flux density distributions are presented in Fig. 11. 
For the range of model parameter values, i.e., Ej/E c = 1CP 2 — 
10 2 , T C /3 C = 0.5 - 20, e B = 10~ 2 - 10" 4 , t F = 10 2 - 10 4 
s and E c = 10 50 erg, the optical peak flux density f Vopt ranges 
from ~ 0.01 /iJy to ~ 0.1 Jy, and the X-ray peak flux density 
is f„ x m 0.001 /iJy —1 mjy for z — 2. For typical parameter 
values, i.e., Ej/E c = 10~\ T c p c = 3, e B = 10~ 2 and t F = 10 2 
s, the fluxes are f Vopt ~ 0.1 mjy, f„ x ~ 1 /iJy. The ratio of 
the peak flux densities at optical and X-rays is roughly constant: 



/"opt / ff 



10 , since the optical band is much closer to the 



spectral peak than the X-rays. The peak flux density is highest for 
r c w 5 — 10, i.e., when the cocoon is mildly relativistic. 

The peak flux densities have very broad ranges. The higher 
ends of the ranges are high compared to the afterglows, but they 
correspond to some certain extreme model parameter value (i.e., 
Ej/E c ~ 10 2 - 10 3 ) and higher values among model parameter 
ranges (e.g., r c ~ 5 — 10). The flux values corresponding to typ- 
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Figure 11. The contours of the observed peak flux densities from the late jet - cocoon interaction. Left panels: at the optical band (v = 2 eV). Right panels: at 
the X-ray band (u = 1 keV). The flux density values labeled on the contours are the demarcation values for two neighbouring contour belts. The ratio of the 
late jet's duration over its delay t^ ur /tp is 0.3 and the redshift is z = 2. Top panels: for tp = 10 2 s; other parameter values are same as in Fig. 3. Middle 
panels: same as in top panels except for eg = 10 -4 . Bottom panels: same as in top panels except for tp = 10 4 s. 



ical model parameter values are comparable to those of observed 
afterglows (see Figs. 17 - 18). 

Fig. 11 also shows the dependence of peak flux densities on 
es and tp. When e_e varies from 10~ 2 to 10 -4 , /„ decreases by 
a factor of ~ 10. This reflects the fact that synchrotron electrons' 
peak specific radiation power and characteristic frequencies are all 
linearly dependent on B' . When tp increases from 10 2 s to 10 4 
s, f v drops by a factor of ~ 10 2 . This is because the interaction 
radius n oc tp and therefore both the cocoon and jet densities are 
smaller at n and so is B' for a larger t p . However for r c /3 c < 2, f v 
increases for a larger tp\ this is due to the optical thickness when 
tp is small (see Sec. 4.5.1). 

Figs. 12 and 13 show the distribution of peak times of the light 



curves, according to which the peak time can be approximated by 
tp, particularly for large values of tp. We measure the full width 
at half maximum (FWHM) of the light curves to characterize the 
pulse width, whose distribution is shown in Figs. 14 and 15. The 
ratio of the pulse width to the peak time, At/t, has a broad range 
of 0.01 - 0.5; typically, for lower T c f3 c or Ej/E c , pulses are wider. 
The optical pulses are often wider than X-ray pulses, but the differ- 
ence is not large. 
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Figure 12. Contours of the light curve peak time for the emission from the 
late jet - cocoon interaction, since the burst trigger and normalized by the 
delay time of the late jet. Top: for the optical light curve; Bottom: for the X- 
ray light curve. Model parameter values are tdur/tF = 0.3 and tp = 100 
s. Other parameter values are same as in Fig. 3. 



5 OBSERVATIONAL IMPLICATION AND DETECTION 
PROSPECTS 

5.1 The Late jet - SNa ejecta interaction 

We find that late jet - SNa ejecta interaction produces a thermal 
transient due to the breakout of a "late" cocoon produced by the late 
jet crossing the SNa ejecta. Thermal X-ray emission was detected 
from X-Ray Flash (XRF) 060218 (Campana et al. 2006). But its 
long lasting, slowly variable light curve suggests that it originated 
in the shock breakout of a quasi-spherical, mildly relativistic ejecta 
and the late jet scenario does not apply to this event. For the prevail- 
ing X-ray flares detected in GRBs, spectral fit shows no compelling 
evidence for a thermal component (Falcone et al. 2007). 

The non-detection of a thermal X-ray transient in GRBs with 
flares indicates that the late jet - SNa ejecta interaction is weak 
or non- existent. There are three possible reasons for the non- 
detection: (i) The cavity in the polar region of the progenitor star 
created by the main GRB jet is still open, and in that case the in- 
teraction between the late jet and the partially filled cavity is weak 
resulting in a significantly lower signal than we have estimated in 
Sec. 3.3. The cavity can be kept open, either because the time was 
too short for the cavity to fill up, or due to a continuous, low-power 
jet that precedes the late jet. This low-power jet might also keep 
the cavity in the cocoon open (see Section 4.2). (ii) The cocoon 
that is in front of the "late" cocoon and the late jet could be opti- 
cal thick, e.g., t c ~ 0.1 E c ^iY~'lt F 2 (also see Figs. 5 - 6), and 
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Figure 13. Same as Fig. 12 except for tp = 10 4 s. 



block the thermal transient signal, (iii) It is also possible that even 
though the GRB was associated with the death of a massive star 
the stellar envelope collapsed producing no supernova at all. In the 
failed SN explosion the whole stellar envelope would collapse on a 
free-fall time scale of a couple of hundred seconds (see, e.g., Ku- 
mar, Narayan & Johnson 2008), form a torus and then be accreted. 
Note that even in this extreme situation, the cocoon associated with 
the main jet should still exist because the cocoon is created when 
the head of the prompt GRB jet passes through the stellar envelope, 
long before the envelope reacts to the core collapse. 



5.2 Late jet - cocoon interaction 

A schematic sketch of light curves due to emission from the late jet 
- cocoon interaction, superposed on the underlying external-shock 
afterglow component, is illustrated in Fig. 16. The expected emis- 
sion from late jet - cocoon interaction has the following features: 

(1) Peak flux density. The peak flux has a broad distribution 
(see Fig. 11). Except for the cases of a very slow cocoon or a very 
low jet-to-cocoon energy ratio, the emission is fairly bright. For 
instance, the range of the calculated X-ray peak flux densities cor- 
responds to a flux in the 0.3-10 keV band of ~ 10" 14 — 10~ 8 erg 
s" 1 cm" 2 (for a spectral index /3x = —1-1) while the Swfit / XRT 
sensitivity is 2 x 10" 14 erg s" 1 cm" 2 in 10 4 s (Gehrels et al. 2004). 
Therefore, the X-ray emission from this interaction is detectable by 
Swift / XRT for most of the parameter space. 

(2) A small At/t (< 0.5; see Figs. 14 and 15). 

(3) A non-thermal spectrum. The emission is mostly non- 
thermal, except when the cocoon speed is sub-relativistic and the 
delay of the late jet is small (e.g., tp <C 10 3 s) thus the thermal 
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Figure 14. Contours of the full width at half maximum (FWHM) of the 
light curve from the late jet - cocoon interaction, normalized by the peak 
time of the light curve since the burst trigger. Top: for the optical light curve; 
Bottom: for the X-ray light curve. Model parameter values are td ur /tp = 
0.3 and tp = 100 s. Other parameter values are same as in Fig. 3. 



photospheric emission might dominate. The X-ray band is in the 
"fast cooling" spectral regime when the emission is non-thermal, 
and the optical is near the synchrotron self absorption frequency. 
The flux density ratio of optical and X-ray is roughly ~ 10 2 for 
p — 2.5 (this ratio is larger for larger values of p). This implies that 
whenever a X-ray pulse is observed, we expect to see an accompa- 
nying increase of the optical flux. 

Could some of the late X-ray flares in GRBs originated from 
the late jet - cocoon interaction? We compare the properties of these 
flares with the prediction of our calculations in the following. 

(1) The observed peak count rate in the XRT band (0.3 - 10 
keV) for flares is distributed over the range of 0. 1 - 100 counts s~ 
(Chincarini et al. 2007). Using the empirical instrument conversion 
factor this translates to f vx ~ l^tJy — lmjy. This range for the 
observed X-ray flare flux is roughly what is expected for the late jet 
- cocoon interaction (cf. Fig. 11). 

(2) The observed value for At/t lies in the range of 0.02 to 
0.6, with a mean value of XRT band (0.3-10 keV) 0. 1 (Chincarini 
et al. 2007). This range is consistent with our calculation (see Figs. 
12 and 13). 

(3) It is rare to find an optical flare accompanying a X-ray 
flare. This is in part due to the fact that very few simultaneous opti- 
cal observations were made in most cases. Nevertheless, we do find 
four cases where an optical flare or rebrightening is reported and the 
X-ray data during the time either is missing or does show a flare. 
These four cases - GRB 050904 (shows an optical flare but its X- 
ray covrage was too sparse to identify a simultaneous flare; Boer et 
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Figure 16. The schematic light curve for the emission from the late jet 
- cocoon interaction superposed on the underlying afterglow light curve. 
Note the zero time in this figure is the burst trigger time, different from the 
zero times in Figs. 8 - 10. 



al. 2006), 060206 (Stanek et al. 2007; Wozniak et al. 2006), 060210 
(Stanek et al. 2007; Curran et al. 2007) and 080129 (Greiner et al. 
2009) - are shown in Fig. 17. The last burst shows an early optical 
flare but without a simultaneous X-ray coverage, and a very late 
(2 x 10 5 s since trigger) rebrightening in both optical and X-rays. 
Among them, the achromatic darings in GRB 050904, 060210 and 
080129 (the very late rebrightening in this burst) are the most likely 
candidates for a late jet - cocoon interaction event. 
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Figure 17. Four GRBs that show a late flare or rebrightening in optical afterglow light curve and, in cases where simultaneous X-ray observations are made, 
a contemporaneous X-ray flare. Top left: GRB 050904; filled circles in red are optical data and filled circles in black are X-ray data (from Boer et al. 2006). 
Top right: GRB 060206, and bottom left: GRB 060210; the open and filled circles in red color (with unnoticeable errors) are the optical data (from Stanek et 
al. 2007). Bottom right: GRB 080129 (from Greiner et al. 2009). 



On the other hand, there are three cases in which simultaneous 
optical observations were available but no optical flare was detected 
at the time of very strong X-ray flare (the X-ray flux increased by 
factors ~ 100 in some of these cases), e.g., GRB 060418, 060607A 
(Molinari et al. 2007) and 060904B (Rykoff et al. 2009) which are 
shown in Fig. 18. This shows that not all X-ray flares are due to the 
late jet - cocoon interaction. However, neither the late jet nor the 
cocoon is ruled out in these three X-ray flares. The lack of an optical 
flare indicates that in these cases it is probable that a low-power, 
continuous jet with a luminosity of Lj,i ow (t) = 2 x lO 48 ^ 2 erg 
s" 1 (total energy ~ a few x 10 49 erg) has kept the cavity open (see 
Section 4.2) and the internal shocks between the more powerful late 
jet and the preceding, slower, low-power jet gave rise to the X-ray 
flares. 



6 SUMMARY 

Observations of X-ray flares in many GRB afterglows suggest the 
existence of a late jet from a long-lived central engine of a GRB at 
~ 10 2 s but possibly 10 4 - 10 5 s after the main GRB event. Adopt- 
ing the collapsing massive star origin for long-duration GRBs, and 
assuming that the supernova explosion to be at approximately the 
same time as the GRB, we have investigated the interactions of this 
late jet with the SNa ejecta and, with a cocoon that was left behind 
when the main GRB jet traversed the progenitor star. 

We find that late jet - SNa ejecta interaction should produce a 
thermal transient, lasting about ~ 10 s and with a peak photon en- 
ergy at a few keV, that should precede or accompany the flare. The 
luminosity of this transient is proportional to the kinetic luminosity 
of the late jet and can be as high as a few x 10 48 erg s" 1 . The lu- 
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minosity is smaller if the polar cavity created by the main GRB jet 
is only partially filled. This thermal transient is similar to the one 
associated with the breakout of the main GRB jet. Although it has 
a lower luminosity, its later occurrence makes it easier to detect. 
The observation of this signal can provide another evidence for the 
massive-star origin of GRBs and new information on the GRB - 
SNa association. 

The fact that no such thermal transient was observed so far 
implies that the late jet - SNa ejecta interaction is suppressed. This 
could happen if the polar region of the ejecta was evacuated by the 
main GRB jet and cavity has not sufficiently refilled itself, espe- 
cially for a not-too-late jet, e.g., ~ 10 2 s, or the cavity could 
be kept open by a continuous, low-power jet. A strong thermal 
transient signal can also be blocked by the optically thick cocoon, 
which would be the case if the cocoon is slowly moving. An al- 
ternative possibility is of a failed supernova scenario in which the 
entire stellar envelope collapses in a free fall time of a few x 10 2 
s. In this case there is no supernova associated with the GRB and 
the late jet, if it is late enough, does not have to cross the stellar 
envelope. 

The late jet interaction with the cocoon would cause a flare or 
rebrightening, superposed on the afterglow light curves, at both the 
optical and X-ray bands. This flare would have a pulse-width-to- 
time ratio At/t < 1 (the expected distribution of At/t is similar 
to that for X-ray flares). Depending on model parameters, we find 
for a burst at a redshift z = 2 that the peak flux density at optical 
fvopt ranges from 0.01 /xjy to 0.1 Jy (V-band apparent magnitude 
29 to 1 1.5) and at X-rays /„ x ranges from 0.001 /xJy to 1 mJy. For 
typical parameters f Vopt ~ 0.1 mJy (V-band magnitude ~ 19) and 
f Vx ~ 1 ^tjy. Observational identification of this emission would 
verify the existence of the cocoon produced when the GRB jet tra- 
versed the progenitor star, thus it would be another confirmation of 
the collapsar model for long duration GRBs (e.g., MacFadyen & 
Woosley 1999; Ramirez-Ruiz et al. 2002; Matzner 2003; Zhang et 
al. 2004). 

The late jet - cocoon interaction might have already been de- 
tected in four GRB afterglows in which simultaneous X-ray and 
optical flares with At/t <C 1 were observed after the prompt emis- 
sion has died off (see Fig. 17). From those candidate events, one 
can learn about the energetics of late jet and the cocoon by utilizing 
the emission calculation presented in this paper. Let us consider the 
flare event in GRB 050904 as an example. We find the most proba- 
ble model parameters - for this burst at 2 = 6.3 and with ti? = 70 
s - that produce the observed peak /„ t and f Vx (data from Boer 
et al. 2006; Cusumano et al. 2007; Gou, Fox & Meszaros 2007) to 
be E c « 10 51 erg, r c « 20 - 50, Ej w 10 52 erg and Vj « 500. 
Those high energetics seem consistent with the very luminous na- 
ture of both the burst and the flare. 

There are three cases in which no optical flare was detected at 
the time of a strong X-ray flare, even though a number of optical 
telescopes have been observing these bursts at the time of the X- 
ray flares (Fig. 18). This shows that not all X-ray flares are due 
to the late jet - cocoon interaction. However, neither a late jet nor 
the cocoon can be ruled out in these cases. It is possible that a 
low-power jet preceding the late jet with a total energy of at least 
10 49 erg had kept the cavity in the cocoon open, so that the late jet 
- cocoon interaction was suppressed. If correct this implies a low 
level continuous emission from the central engine at the level of ~ 
10 47 (t/lOs) -2 erg s -1 lasting for ~ 10 2 s, assuming a radiation 
efficiency of ~ 0.1. 



10 5 
10" 
•3 10 3 

*? 10 s 

w 

I 10' 
x 

E 10° 
10-' 

Id" 8 



1 1 1 mi 1 — 1 — 1— r 




n — 1 — m_ 

GRB 060418 



■ t=130 s 



• H » z' 

♦ J • X 



Lll_ 



10 2 10 3 10" 10 5 

Time since burst (s) 




10 a 10 3 10" 

Time since burst (s) 



10 5 



io- 7 p%^ 

10 ! 



E 

o 



1(T 



g 10- 1 



10" 1 



■. / 



GRB 060904b 
z=0.70 



1 10 100 1000 10000 

Time Since Burst (s) 

Figure 18. The three GRBs that show prominent late X-ray flares but with- 
out simultaneous optical flare apparent in the afterglow light curve. Top: 
GRB 060418; Middle: 060607a (both from Molinari et al. 2007). Bottom: 
GRB 060904b; blue triangles are BAT data extrapolated to X-ray band, ma- 
genta squares are XRT data and red circles are optical data (from Rykoff et 
al. 2009). 
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